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Our  long-term  objective  Is  an  understanding  of  the  molecular  basis  of 
calmodulin  activation  of  target  enzymes.  We  have  chosen  the  Interaction  of 
calmodulin  with  peptides  as  an  Initial  paradigm  to  simplify  the  study  of 
the  molecular  details  of  the  peptide-protein  Interaction.  Our  Initial 
focus  has  been  the  development  of  methodology  for  the  experimental  study  of 


the  molecular  dynamics  of  calmodulin  binding-peptides,  and  of  calmodulin- 


peptide  complexes.  In  particular  an  evaluation  of  the  conjoint  use  of 


fluorescence  and  NMR  spectroscopy  In  such  studies. 


DETAILED  PROGRESS  REPORT 


Most  molecular  dynamics  stimulations  of  protein  structure  have  so  far 


focused  on  events  occurring  In  the  low  picosecond  regime  -  <200  picosecond 
(1).  The  results  of  such  simulations  have  been  Interesting  and  provoca¬ 


tive,  but  they  need  corroboration  from  experimental  data  to  be  validated. 


However,  quantitative  measurement  of  picosecond  motions  is  difficult.  Both 


NMR  and  fluorescence  spectroscopy  are,  in  principle,  able  to  gain  access  to 


this  time  domain  and  to  provide  quantitative  data,  but  to  date  it  has  been 


difficult  to  use  these  two  experimental  methods  to  study  motions  In  preci¬ 


sely  the  same  system.  Conjoint  use  of  both  methods  would  allow  the  data 


from  one  experimental  technique  to  corroborate  (or  refute)  the  Inferences 


drawn  from  the  data  of  the  other. 


Our  basic  experimental  approach  has  thus  been  to  develop  an  appropriate 


paradigm.  To  achieve  this  goal  we  have  synthesized  a  tryptophan  molecule 


with  C  enrichment  In  the  c6i  position  of  the  Indole  ring  (2), 
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Subsequently,  we  have  Incorporated  this  tryptophan  Into  two  sets  of  pep¬ 
tides  -  first  a  set  of  small  hydrophobic  peptides,  and  second  Into  melittin 


and  a  series  of  melittin  analogs.  These  peptides  have  then  been  examined 
by  the  use  of  ^’c-NMR  and  fluorescence  spectroscopy.  The  results  have 
shown  convincingly  the  value  of  this  approach. 

(I)  Details  of  the  Initial  study  with  the  small  hydrophobic  peptides  are 
provided  In  the  manuscript  attached  to  this  report  and  so  won't  be  con¬ 
sidered  here.  A  simple  summary  suffices. 

1.  The  measurements  of  Ti,  T2,  and  NOE  for  the  ”c-H  vector  of  the 
^*C-enr1ched  tryptophan  appears  to  provide  realistic  values  of  the 
rotational  correlation  times  of  the  tryptophan  side  chain.  This 
Inference  Is  based  on  the  similarity  of  the  times  calculated  from  the 
NHR  and  fluorescence  data. 

2.  The  rotational  correlation  times  of  the  Indole  side  chain  of  the  pep¬ 
tides  In  solution  are  generally  In  excess  of  100  ps.  These  values  are 
1  to  2  orders  of  magnitude  different  from  those  predicted  from  molecu¬ 
lar  dynamics  simulations  of  peptides  {In  vacuo)  and  suggest  that 
meaningful  time  constants  from  such  simulations  are  likely  to  be 
obtained  only  If  these  simulations  are  done  with  the  peptides  (or  for 
that  matter  proteins)  In  a  solvent  environment. 

3.  Detailed  mathematical  analysis  of  the  NMR  data  have  thus  far  relied 
exclusively  on  the  use  of  the  Liparl  and  Szabo  model  (3).  The  results 


cc 


have  shown  the  inherent  complexity  of  the  analysis  In  that  the  local 
rotational  (or  llbratlonal)  rate  of  a  molecule  (e.g.,  tryptophan  side 
chain)  calculated  by  use  of  the  Liparl  and  Szabo  formalism  (xe)  Is  very 
sensitive  to  the  order  parameter(S)  and  to  the  overall  rotational 
correlation  time  (I.e.,  Xm)  of  the  substrate.  Consequently,  the 
accuracy  of  the  NHR  measurement  Is  constrained.  Fortunately,  the  error 
Is  unlikely  to  be  even  as  great  as  an  order  of  magnitude,  and  with 
fluorescence  anisotropy  data,  especially  time-dependent  fluorescence 
anisotropy  to  buttress  the  NMR  data,  reasonable  correlation  times  will 
be  realizable.  However,  a  different  formalism  seems  needed  par¬ 
ticularly  for  the  analysis  of  the  dynamics  of  flexible  chains.  We  are 
currently  examining  the  possibility  that  an  appropriate  theory  might 
appear  from  treatment  of  the  peptide  as  a  flexible  polymer  whereupon 
Flory-type  polymer  dynamics  may  be  applied  to  assess  dynamics.  As  we 
noted  above,  the  details  of  the  Initial  study  are  In  the  manuscript 
attached. 

(II)  To  continue  the  assessment  of  ^’C-NMR  as  a  method  for  studying  peptide 
and  protein  dynamics,  we  have  synthesized  and  purified  the  following  pep¬ 
tides. 

*  * 

1.  Helittin  (MLT)  G1GAVLKVLTTGLPALISWIKRKRQQCONH2 

*  * 

11.  Me11tt1ni7  G1GAVLKVLTTGLPALWSILKRKRQQCONH2 

** 

111.  Helittinii  G1GAVLKVLTWGLPALISLIKRKRQQCONH2 

*  * 

1v.  Mel  11  ting  G1GAVLKVWTTGLPALISLIKRKRQQCONH2 


In  each  of  these  peptides  two  ”C-enr1chment  sites  were  Included  -  one  on 
the  a-carbon  of  glycine  at  position  13  In  the  sequence  and  the  other  at  the 
cSi  position  of  the  arotnatlc  ring  (as  described  above).  All  of  the  analogs 
of  melittin  showed  properties  akin  to  that  of  the  parent  peptide.  Most 
notably,  (1)  they  all  showed  a  strong  tendency  to  form  a-hellces  In  hellco- 
genlc  solvents  such  as  methanol  or  trifluoroethanol ;  (11)  they  all  lyse  red 
blood  cells;  (111)  they  all  bind  to  calmodulin.  However,  there  were  also 
some  notable  differences  from  melittin.  The  most  significant  among  them 
being  (1)  the  marked  variation  In  cytolytic  potency  (rank  potency  order 
MLT>MLTi7>MLTii»MLTg) ;  and  (11)  the  failure  of  any  of  the  analogs  to  form 
a  tetramer.  By  use  of  attenuated  total  reflectance  FT-1R  procedures  we 
have  shown,  moreover,  that  the  pattern  of  Interaction  of  the  analog  pep¬ 
tides  with  lipids  varies.  The  apparent  lipid  bilayer  penetrability  shows 
the  same  rank  order  among  the  peptides  as  for  cytolysis.  The  greatest  dif¬ 
ference  Is  seen  for  MLTg  which  appears  to  reside  primarily  on  the  lipid 
surface  and  displays  a  partially  helical /partially  B  sheet  configuration  - 
c.f.,  melittin  whose  amino  terminal  dodecapeptide  appears  to  penetrate  Into 
the  lipid  bilayer  as  a  helix  whose  long  axis  lies  normal  to  the  plane  of 
the  b1 layer  (4). 

The  ^’C-enr1ched  melittin  and  melittin  analogs  therefore  provide  us  with  an 
excellent  model  for  studying  basic  features  of  the  ^’c-NMR  approach  to 
quantitation  of  molecular  motion.  The  *’c-a  carbon  group  of  glycine,  kept 
In  a  comnon  location  In  all  of  the  peptides,  provides  direct  information  on 
the  dynamics  of  the  peptide  main  chain  and,  therefore,  allows  assessment  of 
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the  overall  dynamics  of  the  peptide.  The  *’C-enr1ched  Indole  ring  provides 
data  almost  exclusively  on  the  mobility  of  the  aromatic  side  chain.  We 
have  now  studied  In  detail  the  solution  dynamics  of  all  four  peptides  under 
conditions  where  we  expect  the  peptides  to  be  (1)  In  random  (coll)  con¬ 
figuration;  (11)  existing  as  monomeric  helices  (e.g..  in  70%  methanol /water 
mixtures,  V/V);  and  (111)  tetrameric  -  I.e.,  In  solutions  with  high  salt 
concentration.  The  measurements  have  been  very  productive.  So  far,  data 
have  been  analyzed,  as  described  above,  by  the  use  of  the  Liparl  and  Szabo 
"model -free*  approach  as  Implemented  by  Weaver  et  a1.  (see  manuscript 
attached).  The  results  and  Inferences  are  briefly  sunnarlzed  below. 

(I)  For  the  random  coll  peptides,  the  apparent  rotational  correlation 
times  calculated  from  T^,  and  NOE  measurements  show  an  expected, 
distinct  difference  between  the  motions  of  the  ”c-g1yc1ne  site  and  the 
**C-tryptophan  (detailed  values  are  provided  In  Table  1).  The  motion 
of  the  Indole  side  chain  showed  rotational  correlation  times  of 

ca  ZOO  ps.  for  the  peptides  dissolved  In  0MSO-d5:02O  (90:10  V/V) (see 
table  I),  whereas  the  glycine  a-carbon  showed  correlation  times  >  5 
nsec  even  for  the  peptides  In  free  solution  (as  monomers?) 

(II) An  Increase  In  salt  concentration  Is  known  to  make  HLT  helical  and  to 
provoke  formation  of  a  melittin  tetramer.  The  ”C-NMR  results  showed 
clearly  the  evolution  of  this  tetrameric  structure  as  the  Ionic 
strength  was  Increased  In  a  solution  of  the  ^’c-enr1ched  melittin.  The 
process  was  reflected  In  a  progressive  up-fleld  shift  in  the  position 


of  the  ^’C-resonance  of  both  the  tryptophan  side  chain  and  the  glycine 
a-carbon. 

Interestingly,  the  NMR  data  also  suggested  the  presence  of  oligomeric 
species  other  than  tetramer  (most  likely  melittin  dimers).  The  formation 
of  the  tetramer  was  also  reflected  In  a  tenfold  Increase  In  the  apparent 
rotational  correlation  time  of  tryptophan  side  chain  motion  and  a  similar 
change  In  correlation  time  for  the  glycine  a-carbon.  However,  Increased 
Ionic  strength  failed  to  significantly  affect  the  ”c-NMR  spectra,  Tj,  or 
NOE  of  the  ^’c-resonances  In  any  of  the  melittin  analogs  supporting  the 
view  that  these  analogs  do  not  form  oligomers. 

The  changes  In  chemical  shift  for  both  the  tryptophan  and  glycine 
^’c-resonances  support  the  Inference  that  both  carbon  moieties  move  into 
less  polar  environments  as  melittin  becomes  a  tetramer.  This  Inference  Is 
supported  by  examination  of  the  x-ray  crystal  structure  of  the  tetramer. 
The  Increases  In  rotational  correlation  time  stem  from  the  effective 
Increase  In  molecular  volume  attending  oligomer  formation,  and  from  the 
Increased  restriction  on  motion,  as  both  the  tryptophan  side  chain  and  the 
glycine  residues  are  part  of  the  peptide-peptide  Interaction  surfaces  in 
the  tetramer.  The  ^’c-NMR  results  on  tryptophan  side  chain  behaviour  are 
strongly  corroborated  by  the  results  of  fluorescence  measurements.  A 
manuscript  describing  this  work  Is  currently  being  prepared  for  publica- 


The  data  we  have  gathered  during  the  past  year  have  thus  established 
the  usefulness  of  the  combined  ^*C-NMR  and  fluorescence  approach  to  study 
of  peptide  dynamics.  The  process  needs  now  to  be  applied  to  more  complex 
systems  (see  plans  for  grant  year  1987-1988). 

(Ill)  Molecular  Dynamics  Simulations 

As  stated  earlier,  one  of  the  principal  long-term  objectives  of  this 
research  Is  the  study  of  the  Influence  of  protein  dynamics  on  the  function 
of  calmodulin  as  a  Ca*'~''-dependent  regulator  of  enzyme  function.  Work  on 
this  aspect  of  the  project  has  only  just  now  begun  because  of  the  lack  of  a 
crystal  structure  for  calmodulin  I.e.,  lack  of  the  x-ray  crystal  coor¬ 
dinates.  (The  crystal  structure  of  calmodulin  Is  known,  but  the  coor¬ 
dinates  have  yet  to  be  deposited  In  the  Brookhaven  data  bank.)  During  the 
past  two  months  we  have  attempted  to  surmount  this  problem  In  the  following 
way.  The  a-carbon  backbone  (main  chain)  coordinates  are  available.  Since 
the  gross  secondary  form  of  the  protein  Is  known.  It  should  be  possible  to 
use  the  molecular  dynamics  algorithm  with  energy  minimization  procedures  to 
reconstruct  the  side  chain  orientations.  The  problems  Inherent  in  such 
reconstructions  are,  of  course,  well  known  to  those  who  have  tried  to  fold  a 
protein  of  known  amino  acid  sequence  Into  a  secondary/tertlary  structure. 

In  our  case  the  difficulty  of  the  latter  procedure  Is  lessened  considerably 
since  the  secondary  structure  Is  essentially  already  provided  by  the  con¬ 
figuration  of  the  a-carbon  chain.  However,  to  give  the  greatest  credence 
to  the  results  of  any  calculations  we  might  do  on  calmodulin,  we  have  first 


secured  the  complete  coordinate  set  for  the  structurally  homologous 
troponin-C  (turkey).  These  have  been  kindly  supplied  by  Dr.  M.  N.  James  of 
the  University  of  Alberta.  We  have  reduced  the  data  on  troponin-C  to  the 
a-carbon  set  and  are  currently  using  energy  minimization  and  molecular 
oynamlcs  simulations  with  CHARMM  to  "build”  the  troponin-C.  The  calculated 
structure  will  then  be  compared  with  that  from  the  x-ray  coordinates  (see 
"experimental  plans  below). 

EXPERIMENTAL  PLANS  FOR  1987-1988 
A.  NMR  and  Fluorescence  Measurements 

1.  Now  that  the  studies  on  peptides  In  solution  have  been  completed, 
we  will  begin  analysis  of  the  calmodulin-peptide  complexes.  We 
know  that  melittin  and  the  three  analogs  described  above  all  bind 
to  calmodulin.  Also,  we  know  from  other  studies  that  these  pep¬ 
tides  assume  a  (presumably  amphiphilic)  helical  form  when  bound. 

Our  objectives  with  the  “C-NMR  measurements  will  be  to  determine 
the  dynamics  of  the  a-carbon  of  the  ”C-glyclne  and  of  the  Indole 
ring  In  the  complexes,  under  conditions  In  which  (1)  the  tryptophan 
forms  part  of  the  Interacting  surface  (I.e.,  lies  on  the  apolar 
face  of  the  amphiphilic  helix)  and  should,  therefore,  be  tightly 
"packed"  In  the  overall  structure  of  the  protein,  peptide  complex, 
or  (11)  the  tryptophan  Is  relatively  free  to  rotate  I.e.,  when  the 
Indole  side  chain  lies  on  the  hydrophilic  surface  of  the  CAM-bound 
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peptide  (e.g.,  MLTn).  In  all  cases,  the  ”C-glyc1ne  resonance 
should  show  highly  restricted  motion  in  the  complex  i.e.,  should 
display  an  effective  order  parameter  of  1. 

2.  It  will  be  necessary,  concomitantly  with  the  NMR  measurements,  to 
re-examine  the  reliability  of  the  Lipari  and  Szabo  model  for  calcu¬ 
lation  of  dynamic  parameters.  Our  studies  to  date  have  shown 
clearly  the  limitation  of  the  model  under  conditions  where  the 
overall  rotational  correlation  time  (1^  in  our  notation)  cannot  be 
specified  i.e.,  has  not  been  corroborated  by  data  from  some  other 
independent  measurement.  The  problem  is  that  tm  is  not  easily 
measured.  However,  in  a  folded  protein  whose  molecular  dimensions 
are  known,  a  reasonable  estimate  of  can  be  calculated  by  use  of 
the  Stokes-Einstein  equation.  Moreover,  time-dependent  anisotropy 
data  should  provide  another  estimate  of  Xm  which  one  hopes  would  be 
similar  to  that  calculated  from  the  molecular  dimensions.  Thus, 
all  NMR  measurements  will  continue  to  be  complemented  with  deter¬ 
minations  of  time-dependent  fluorescence  anisotropy  decays.  For 
the  latter  our  Nd:YAG-pumped,  picosecond  resolving  multi  frequency 
phase  fluorometer  will  be  used.  These  measurements  will  be 
employed  for  the  assessment  of  the  proteins'  rotational  correlation 
times.  Unfortunately,  studies  of  peptide  dynamics  do  not  enjoy 
such  an  advantage.  A  detailed  study  of  polymer  dynamics  will  be 
undertaken  during  the  coming  year  to  determine  possible  applicabi¬ 
lity  to  the  analysis  of  the  NMR  data. 
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B.  Molecular  Dynamics  Calculations 

In  work  supported  by  other  grants,  we  have  completed  molecular  dynamics 
studies  of  cyclosporine,  scorpion  neurotoxin  variant  3  and  ribonuclease 
Ti*  The  program  CHARMM  was  used.  The  experience  gained  has  con- 
sidp-ably  facilitated  our  Initial  attempts  to  “build"  TnC  from  a-carbon 
coordinates.  We  expect  to  have  solved  the  problems  with  this  construc¬ 
tion  within  a  few  weeks  of  the  new  grant  year.  The  next  stage  Is  to 
build  calmodulin  from  Its  a-carbon  coordinates  and  to  proceed  with 
molecular  dynamics  calculations.  The  latter  will  not  be  straightfor¬ 
ward  for  several  reasons.  First,  the  x-ray  structure  was  done  on  the 
molecule  at  acidic  pH.  All  available  evidence  suggests  that  under 
such  conditions,  CAM  has  a  significantly  different  configuration  from 
that  existing  at  neutral  pH.  Moreover,  the  protein's  Ca++  binding  pro¬ 
perties  perforce  are  different  at  acidic  pH.  The  second  problem  Is 
strongly  related  to  the  first,  namely  that  the  high  anionic  charge 
complicates  our  modeling  because  It  requires  that  we  account  carefully 
for  a  high  electrostatic  charge.  In  our  view,  this  negates  a  simple 
“vacuum"  simulation  and  necessitates  inclusion  of  water.  Even  with  the 
latter,  the  electrostatic  problem  Is  ameliorated  but,  we  realize,  cer¬ 
tainly  not  solved.  Thus,  we  know  from  CD  measurements  that  the  confor¬ 
mational  response  of  CAM  to  the  binding  of  Ca++  Ions  Is  different  If 
the  protein  Is  dissolved  In  media  of  low  or  "physiological"  Ionic 
strengths.  We  cannot  currently  correct  the  simulations  for  Ionic 
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effects  but  would  argue  that  simulations  in  an  aqueous  environment 
should  provide  at  least  a  good,  if  preliminary,  assessment.  We  will 
most  likely  employ  a  periodic  boundary  calculation  in  order  to  include 
water. 


Once  we  are  satisfied  with  the  results  of  our  construction  of  the  CAM 
structure  we  will  proceed  with  i)  a  calculation  of  a  likely  structure 
for  CAM  in  its  fully  Ca++- replete  form;  ii)  a  molecular  dynamics  simu¬ 
lation  of  the  CAM;  iii)  a  molecular  graphics/molecular  dynamics 
assessment  of  the  likely  peptide  binding  site  on  the  Ca++  replete  mole¬ 
cule.  Clearly  item  number  ii  will  present  the  most  formidable 
challenge  and  the  fate  of  item  iii  is  utterly  dependent  on  the  success 
of  the  initial  dynamics  studies.  CAM  is  large  enough  that  a  simulation 
of  the  entire  molecule  for  a  sufficiently  long  time  and  inclusive  of 
water  is  probably  impracticable.  An  initial  approach  will  be  to  simu¬ 
late  each  major  domain  separately  and  then  try  to  reconstruct  the 
entire  protein. 


The  cost  in  computer  and  personnel  time  will  be  high.  Dr.  Chris 
Haydock,  who  has  become  expert  in  the  use  of  molecular  dynamics  calcu¬ 
lations  will  be  fully  employed  on  the  project.  Simulations  will  con¬ 
tinue  to  be  performed  on  CRAY  computers. 


